MATHEMATICAL ENGINEERING TECHNICAL REPORTS Calculation of Orthant Probabilities by the Holonomic Gradient Method
نویسندگان
چکیده
We apply the holonomic gradient method (HGM) introduced by [9] to the calculation of orthant probabilities of multivariate normal distribution. The holonomic gradient method applied to orthant probabilities is found to be a variant of Plackett’s recurrence relation ([14]). However a straightforward implementation of the method yields recurrence relations more suitable for numerical computation than Plackett’s recurrence relation. We derive some theoretical results on the holonomic system for the orthant probabilities. These results show that multivariate normal orthant probabilities possess some remarkable properties from the viewpoint of holonomic systems. Finally we show that numerical performance of our method is comparable or superior compared to existing methods.
منابع مشابه
Calculation of orthant probabilities by the holonomic gradient method
We apply the holonomic gradient method (HGM) introduced by [9] to the calculation of orthant probabilities of multivariate normal distribution. The holonomic gradient method applied to orthant probabilities is found to be a variant of Plackett’s recurrence relation ([14]). However an implementation of the method yields recurrence relations more suitable for numerical computation than Plackett’s...
متن کاملMATHEMATICAL ENGINEERING TECHNICAL REPORTS Estimation of Exponential-Polynomial Distribution by Holonomic Gradient Descent
We study holonomic gradient decent for maximum likelihood estimation of exponentialpolynomial distribution, whose density is the exponential function of a polynomial in the random variable. We first consider the case that the support of the distribution is the set of positive reals. We show that the maximum likelihood estimate (MLE) can be easily computed by the holonomic gradient descent, even...
متن کاملMATHEMATICAL ENGINEERING TECHNICAL REPORTS Holonomic Gradient Method for the Distribution Function of the Largest Root of a Wishart Matrix
We apply the holonomic gradient method introduced by Nakayama et al. [19] to the evaluation of the exact distribution function of the largest root of a Wishart matrix, which involves the hypergeometric function 1F1 of a matrix argument. Compared to methods based on infinite series expansion in terms of zonal polynomials, whose convergence is slow, the proposed method is more practical for compu...
متن کاملMATHEMATICAL ENGINEERING TECHNICAL REPORTS Holonomic Gradient Method for Distribution Function of a Weighted Sum of Noncentral Chi-square Random Variables
We apply the holonomic gradient method to compute the distribution function of a weighted sum of independent noncentral chi-square random variables. It is the distribution function of the squared length of a multivariate normal random vector. We treat this distribution as an integral of the normalizing constant of the Fisher-Bingham distribution on the unit sphere and make use of the partial di...
متن کاملOptimization of Longwall Panel Location with Regards to the Gradient of Coal Seams (TECHNICAL NOTE)
The paper begins by outlining the role and importance of coal as a source of energy and in the steel industry. It briefly describes the longwall method of working together with the conventional machinery used in the method. A mathematical model is then proposed that shows the relationship between the gradient of the coal seam, that of the face and the entries to the panel. Determination of a mo...
متن کامل